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Abstract. We compare the accurate empirical mass-luminosity (M-L) relation based on five Hyades binary sys- 
tems to predictions of stellar models calculated with various input parameters (helium, metallicity and age) or 
physics (mixing-length ratio, model atmosphere, equation of state and microscopic diffusion). Models based on 
a helium content y~0.28 inferred from the AY/AZ enrichment law are more than 3a beyond the observations, 
suggesting that the Hyades initial helium abundance is lower th an ex pected from its supersolar metallicity. With 



the photometric metallicity ([Fe/H]=0.144±0.013 dex, Grenon 2000) we derive F=0.255±0.009. Because of the 



(y, [Fe/H]) degeneracy in the M-L plane, the uncertainty grows to AY = 0.013 if the metallicity from spec- 



troscopy is adopted ([Fe/H]=0.14±0.05 dex, Cayrel de Strobel et al. 1997). We use these results to discuss the 
Hertzsprung-Russell (HR) diagram of the Hyades, in the (Mv, B — V) plane, based on the very precise Hipparcos 
dynamical parallaxes. Present models fit the tight observed sequence very well except at low temperatures. We 
show that the HR diagram does not bring further constraints on the helium abundance or metallicity of the 
cluster. In the low mass region of the HR diagram sensitive to the mixing-length parameter (omlt), the slope of 
the main sequence (MS) suggests that qmlt could decrease from a solar (or even supersolar) value at higher mass 
to subsolar values at low mass, which is also supported by the modeling of the vB22 M-L relation. We find that 
the discrepancy at low temperatures {{B — V) > 1.2) remains, even if an improved equation of state or better 
model atmospheres are used. Finally, we discuss the positions of the stars at turn-off in the light of their observed 
rotation rates and we deduce that the maximum age of the Hyades predicted by the present models is ~650 Myr. 
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1. Introduction 

Open clusters provide information and strong constraints 
for the stellar evolution theory. They give the opportunity 
to study large numbers of stars spanning a broad range of 
masses and evolutionary stages and that can be assumed 
to have similar age and chemical composition. Depending 
on the cluster studied, observations may give 1) the po- 
sition of the cluster sequence in the HR diagram, 2) the 
density of stars along that sequence, 3) the M-L relation 
if the masses of some binary stars are accessible to ob- 
servation. The analysis of the observational features of a 
given cluster by means of internal structure models allows 
estimating characteristics not directly accessible through 
observation, such as the age or helium content of the mem- 
bers. Furthermore, if the observational data are accurate 
enough, constraints on the physical processes at work in 



the stellar interiors, for instance the various transport pro- 
cesses, can be inferred. 

The Hyades is the nearest moderately rich star cluster. 
It has served for a long time to define absolute magnitude 
calibrations and, in turn, to fix the zero-point of the galac- 
tic and extragalactic distance scales. Also, as an open clus- 
ter, the precise knowledge of its chemical composition and 
age is fundamental for studies of the kinematic and chem- 
ical evolution of our Galaxy. The metallicity [Fe/H] (log- 
arithm of the number abundances of iron to hydrogen rel- 
ative to the solar value) is accessible through photometric 
or spectroscopic observations (see Sect. pj| ). On the other 
hand, the Hyades dwarfs are too cool for helium lines to 
be visible in their spectra and their helium content to be 
determined directly. The helium abundance and age have 
to be derived from the analysis of the observations using 
stellar models. 



Send offprint requests to: Y. Lebreton, e-mail: 
Yveline . LebretonSobspm . f r 



Recently, high-quality observations of the Hyades 
stars have been obtained. Hipparcos data combined with 
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ground-based photometric or spectroscopic observations 
provided a more precise extended HR diagram of the clus- 
ter (Ferryman et al. 1998, Dravins et al. 1997, de Bruijne 
et al. 2000). On the other hand, observations of several 
binary systems in the Hyades yielded a much improved 
M-L relatio n (To rres et al. |1997a| 
& Solensky |l98|; Sderhjelm"^99|) " 

Ferryman et al. ( |1998D, L ebreton et al. ( |1997| ) and 



0.015 



1997b, 1997c; Peterson 



Cayrel de Strobel et al. (1997) analysed the observational 
HR diagram of the Hyades with stellar models and ob- 
tained new estimates of the age, initial helium content 
and metallicity of the cluster (see Sect. 2.1). Recently, de 
Bruijne et al. ( ^000 ) and Castellani et al. ( 2000D compared 
the observational HR diagram with models and discussed 
the different uncertainties, in particular the problems re- 
lated to the color-magnitude calibrations. 

In this paper, we focus on the determination of the 
helium abundance of the Hyades through the analysis of 
both the M-L relation and HR diagram. We show that 
the excellent accuracies on the masses and luminosities 
reached for the Hyades binaries provide rather severe con- 
straints on the helium content of the cluster and allow 
confirming and refining of the helium value previously in- 
ferred from the analysis of the HR diagram. In a second 
step, we discuss our ability to reproduce the various fea- 
tures of the HR diagram cluster sequence with different 
model input physics. 

Sect. H is a brief review of previous studies of the ob- 
served HR diagram of the Hyades. Sect. |^ presents the 
observed Hyades M-L relation. Numerical stellar models 
are presented in Sect. ^. In Sect. ^, we analyse the M-L 
relation by means of stellar models and discuss the im- 
plications for the chemical composition of the cluster. In 
Sect. |[ we analyse and discuss the HR diagram of the 
cluster on the basis of the constraints provided by the M- 
L relation. Summary and conclusions are given in Sect. 

2. HR diagram data and previous analyses 

2.1. Observational data 

Hipparcos measured the positions, proper motions and 
trigonometric parallaxes of ~300 candidate Hyades mem- 
bers allowing new, more detailed studies of the cluster to 
be undertaken. 



Ferryman et al. (1998) assigned the membership of 
218 cluster stars and from the individual trigonometric 
parallaxes of the 134 stars located within 10 pc of the 
cluster center, obtained a distance modu lus (m — M) = 
3.33 ± 0.01 mag. Narayanan & Gould (|l999|) used the 
Hipparcos proper motions to derive statistical parallaxes 
for 43 cluster stars leading to (m — M) = 3.34 ± 0.02 mag. 



Dra vins et al. (1997), and more recently de Bruijne 
et al. ( 2000 ), derived the dynamical parallaxes of ~200 
Hyades members from the relation between the cluster 
space motion, the positions and the projected proper mo- 
tions. As shown in Fig. |l], the two groups, who used the 
same data and similar methods, obtain parallaxes in excel- 
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Fig. 1. Comparison of dynamical parallaxes of Hyades' 
stars as determined by de Bruijne et al. ( 2000D and by 
Dravins et al. ( |1997D : aU stars but one (HIF 28356) have 
their parallaxes within less than 0.4 percent 



lent agreement; the dynamical parallaxes of the stars they 
have in common agree to better than 0.4 percent (except 
for one star, HIF 28356, where the parallax difference is 
about 1 percent). The mean accuracy on the dynamical 
parallaxes is ^^0.5 mas, that is ^3 times better than the 
trigonometric parallax accuracy. 

Fhotometric information can be found in the 
Hipparcos catalogue. The Hipparcos _ffp-magnitude of 
each observed star is given with an accuracy of 0.0015 
mag. The ^-magnitude and [B — V) index of each star 
coming from ground-based measurements are also given 
in the catalogue, the typical accuracy on V is better than 
0.01 mag. 

The mean absolute magnitude accuracy of the Hyades 
stars based on the dynamical parallaxes of Dravins et al. 
( |1997D or of de B ruijne et al. ( |2000D is ctmv ^ 0.05 mag. 
De Bruijne et al. (2000) examined 218 Hyades candidates 
among which they selected 90 secure Hyades members. 
The resulting (My, B — V) colour-absolute magnitude di- 
agram shows a very well-defined MS (i.e. {umv) — 0-05 
mag and {a(^B-v)) — 0.01 mag, see Fig. ^). 

The analysis of high-resolution, high signal-to-noise 
spectra of several Hyades stars has provided precise deter- 
minations of their effective temperatures, Tcs, a nd metal- 
licities, [Fe/H]. Cayrel de Strobel et al. (|l997|) selected 
40 Hyades dwarfs with T^e accurate to typically 50-70 
K; their mean metalhcity [Fe/H]= -1-0.14 ± 0.05 dex, is 



in good agreement with the photometric value [Fe/H] = 
-1-0.144 ± 0.013 dex derived by Grenon (|2000|) from large 
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sets of fiomogeneous observations in the Geneva photo- 
metric system. 

2.2. Previous analyses by means of stellar models 



Cayrel de Strobel et al. ( 1997|) obtained the lower part 
of the cluster MS, in the (Mboi, T'cff) plane, from the po- 
sitions of 40 dwarfs, by combining Hipparcos distances, 
spectroscopic T^s, ^-magnitudes from the Hipparcos in- 
put catalogue and bolometric corrections by Bessel et al. 
(1998). Lebreton et al. (1997) and Ferryman et al. 
(1998) compared the bottom of this sequence, correspond- 
ing to the non-evolved stars, with theoretical zero-age 
main sequences (ZAMS) computed with a Z/X ratio of 
0.034 corresponding to the mean observed [Fe/H]-value 
1^ The ZAMS fitting yielded the initial helium content 
Y= 0.26±0.02 (in mass fraction) and metallicity Z~ 
0.024±0.004 of the cluster, the (quite large) error on 
Y being dominated by the error on the mean value of 
[Fe/H] used ([Fe/H]=-|-0.14±0.05 dexfrom spectroscopy). 

The position of the whole cluster sequence i n the {My, 
B — V) plane was obtained by Ferryman et al. ( 1998D from 
the combination of Hipparcos distances and ground-based 
V"-magnitudes and {B — V) color indexes. The fitting of 
the sequence by model isochrones corresponding to the 
chemical composition of the cluster {Y= 0.26, Z= 0.024) 
yielded an estimate of the age, t ~ 625 Myr. 

In order to obtain a MS as fine as possible, the known 
or suspected unresolved binaries were removed from the 
samples mentioned above and the masses of the few re- 
solved binaries were not used as constraints. Binaries with 
well-determined masses provide valuable constraints for 
stellar models, in particular for the stellar abundances 



through the M-L relation (see e.g. Andersen 1991). Some 
difficulties in fitting parts of the empirical Hyades M-L re- 
lation with models have been encountered by Torres et al. 
( |l997b[ ) and Lastennet et al. ( 1999] ), but the models (avail- 
able from the literature) did not correspond to the cluster 
chemical composition. We now turn to the sample of bi- 
naries well-observed in the Hyades. 



3. Updated data for five binary systems 

We focus on five Hyades binary systems with well- 
determined masses. The data are listed in Table |. The 
dynamical parallaxes, TTdyn, are more precise than the di- 
rect Hipparcos trigonometrical parallaxes, TTtrig (Ferryman 
et al. 199^ ). Therefore, in the following, we have used the 
dynamical parallaxes of de Bruijne et al. ( ^000 ) available 
for the five systems. We also list in Table ^ the TTdyn-values 
derived by Dravins et al. (1997) in their improved version 
(Lindegren 1999, private commimication) . 



^ Z and X are respectively the abundances in mass of 
heavy elements and hydrogen related to [Fe/H] through 
[Fe/H]=log(Z/X) - Ior(Z/X )p, w here {Z/X)q = 0.0245 ac- 
cording to Grevesse & Noels (199J). 



For three systems (51 Tauri, Finsen 342 and 9^ 
Tau), which are double- lined spectroscopic binaries as 
well as visual binaries resolved by speckle, Torres et al. 
(1997a, 1997b, 1997c) derived a complete astrometric- 



spectroscopic orbital solution, and therefore could obtain 
the individual masses and the orbital parallax, TToib, listed 
in Table ^ 51 Tauri is also one of the 25 binary systems 
for which Sderhjelm (199£) could derive the individual 
masses combining Hipparcos data with ground-based ob- 
servations. The resulting masses. Ma — 1.72±0.27 Mq for 
the primary and Mb = 1.31±0.21 M© for the secondary 
are in good agreement with those obtained by Torres et 
al. Also, Martin et al. ( |l998|) obtained the masses of the 
two components of 51 Tau from the analysis of Hipparcos 
data but with a lower accuracy (Ma = 1.756 ± 0.343 M© 
and Mb = 0.953 ± 0.247 Mq). 

For the 6^ Tau system, a single-lined spectroscopic bi- 
nary with astrometric information from speckle and lunar 
occultation measurements, Torres et al. (1997c) obtained 
a partial astrometric-spectroscopic solution where the in- 
formation coming from the velocity amplitude of the B- 
component was lacking. However, it is possible to estimate 
the individual masses, using the parallax obtained inde- 
pendently from Hipparcos data. We give in Table |l| the 
values of the masses of the components of 9^ derived from 
formulae (2), in Torres et al (1997c), using their orbital 
parameters together with the Hipparcos TTdyn. 2 dynamical 
parallax. It is worth noting that 9^ Tau B and (/)342B have 
very similar masses and magnitudes. 

The fifth system, vB22, is a double-li ned sp ectroscopic 
eclipsing binary discovered by McClure ( 1982 ) . Indiv idual 
masses were derived by Feterson & Solensky ( 1988 ) and 
the radii were determined by Schiller & Milone ( 1987 ) . 

The ^-magnitude of each system, the magnitude dif- 
ference between the two components Ay and the corre- 
sponding individual absolute magnitudes My, a and My b 
(corresponding to the TTdyn ,2-values) are listed in Table | 
The individual values of the {B — V) color index of the two 
components of vB22 have been determined by Schiller & 



Milone (1987). For the other systems we infer the individ- 
ual values of {B — V) from the global value of the system 
given in WEBDA (http://obswww.unige.ch/webda/, 
see Mermilliod p_998 ) and from the difference A{B — V). 
For 9^ Tau, A{B — V) has been measured by Feterson 
([1993D . For 51 Tau and 0342, we follow Torres et al. 

V) from the value of AT^ and 
V) relation of 
for the Hyades main sequence 



et al. 



(1997a) who derived A{B - 



of the local slo pe of the empirical {V, B 
Schwan et al 



(|1991|) 



, ^ , For 

the 9^ Tau system, which is composed of a giant and of a 
dwarf, we only know the system (B — T^)-value. 

4. Stellar models and related input physics 



We have calculated stellar models in the mass range 
M e [0.5,3.0] Mq, using the CESAM code (Morel |1997D 
and the input physics described below. 
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Table 1. Parallaxes and masses for the 5 binary systems. The dynamical parallaxes 7rdyn,i are from Lindegren (1999, 
private communication but see also Dravins et al. 1997). The dynamical parallax es TTdy n 2 are from de Bruijne et al. 
(200c ). The Hipparcos trigonometrical parallaxes TTtrig come from Ferryman et al. (1998). As discussed in the text, the 
orbital parallaxes TTorb and individual masses were obtained by Peterson & Solensky ( 198^ ) for vB22 and by Torres 
et al. (1997a, 19971:, 1997c) for the other systems. The radii of the vB22 components are from Schiller & Milone (1987) 
where the probable error has been translated into a mean error. 



name 


HIP 


HD 


■""dyn,! 


'rdyn,2 


TTtrig 


TTorb 


Ma/Mq 


Ra/Rq 








(mas) 


(mas) 


(mas) 


(mas) 


Mb/Mq 


Rb/R© 


vB22 


20019 


27130 


21.15±0.33 


21.16±0.38 


21.40±1.24 


19.8 ±0.4 


1.072 ±0.010 
0.769 ± 0.005 


0.905 ± 0.029 
0.773 ±0.015 


51 Tan 


20087 


27176 


18.29±0.24 


18.31±0.69 


18.25±0.82 


17.9 ±0.6 


1.80 ±0.13 
1.46 ±0.18 




(j> 342 


20661 


27991 


21.27±0.33 


21.29±0.37 


21.47±0.97 


21.44±0.67 


1.363 ± 0.073 
1.253 ± 0.075 




Tau 


20885 


28307 




21.29±0.37 


20.66±0.85 




2.77 ±0.50 
1.28 ±0.13 




Tau 


20894 


28319 


22.21±0.35 


22.24±0.36 


21.89±0.83 


21.22±0.76 


2.42 ± 0.30 
2.11 ± 0.17 





Table 2. Magnitudes and colour indices for the five binary systems. 



name 


HIP 


HD 


V 


Mv,A 


B-V 


{B - V)a 








AV 


Mv,B 


A{B - V) 


(B - V)b 


vB22 


20019 


27130 


8.319 ±0.009 


5.07 ±0.04 




0.713 ±0.017 








2.3 ±0.05 


7.37 ± 0.06 




1.19 ±0.08 


51 Tau 


20087 


27176 


5.65 ±0.01 


2.18 ±0.09 


0.28 ±0.01 


0.25 ±0.05 








1.61 ±0.10 


3.79 ±0.12 


0.19 ±0.01 


0.44 ± 0.02 


(j> 342 


20661 


27991 


6.46 ±0.01 


3.70 ± 0.05 


0.49 ±0.01 


0.466 ± 0.056 








0.34 ± 0.05 


4.03 ± 0.05 


0.05 ±0.01 


0.530 ± 0.020 


Tau 


20885 


28307 


3.84 ±0.01 


0.52 ±0.04 


0.95 ±0.01 










3.50 ±0.05 


4.02 ± 0.07 






Tau 


20894 


28319 


3.40 ± 0.02 


0.48 ± 0.04 


0.18 ±0.01 


0.18 ±0.02 








1.10 ±0.01 


1.58 ±0.05 


-0.006 ±0.005 


0.17 ±0.02 



Equat ion o f state (EOS): the CEFFEOS (Eggleton 
et al. 1973| , Christensen-Dalsgaard |1991| ) has mostly 
been used; it includes the Debye-Hiickel cor rectio ns for 
pressure. The OPAL EOS (Rogers et al. |l996|) that 
includes collective effects has been used at low mass 
for the purpose of comparison. 



Opacities: we used OPAL (Iglesias & Ro gers |1996[) 
complemented by Alexander & Ferguson's (1994) data 
f or T < 10'* K, both being for Grevesse & Noels's 
( |l993| ) solar mixture. 

Therm onuclear reaction rates: Caughlan & Fowler 

Standard mixing-length theory of convection (Bohm- 
Vitense |1958 ). The calibration of the solar model in 
luminosity and radius with a particular set of input 
physics yields the solar mixing-length ratio q;mlt.0 = 



of the nearby stars (Lebreton et al 
Hyades (Ferryman et al. 1998| ) 



1999D and of the 
In addition, the cali- 



bration in luminosity and effective temperature of sev- 
eral nearby visual binaries with determined masses 
yields values of q;mlt close to q^mlt,© (Fernandes 
et al. 1998). From numerical 2-D simulations of con- 
vection in a wide range of T cs an d log g (g is the sur- 
face gravity), Ludwig et al. ( 1999 ) obtained equivalent 
QfMLT-values in the range 1.1-1.8 for solar or subsolar 
[Fe/H]- they did not investigate supersolar metallicites 
as found in the Hyades. For low mass MS stars these 
simulations suggest possible variations of aMLT with 
Toff and log g in the range omlt,© ±0.15. We there- 
fore adopted the solar mixing-length ratio for standard 
models but we discuss the effects of a change of umlt 
as high as ±0.40. 



I/Hp (Hp is the pressure scale- height). The value of 5. Convective core overshooting. We take dov = 0.2Hp as 



omlt.o depends on the external boundary conditions 
adopted. We find omlt, = 1-61 for the Eddington's 
grey model atmosphere and q;mlt, 0=1-79 for Kurucz's 
( |1991D ATLAS9 models (see item 7). Models calcu- 
lated with aMLT,0 also give a good fit of the MS slope 



overshooting distance, a value derived by Schaller et al. 



(1992) from the study of the empirical MS width of 
open clusters. We also discuss the case with no over- 
shooting. 
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6. Microscopic diffusion transports helium and heavy ele- 
ments toward the center of stars and may change their 
surface and core compositions during evolution. It can 
be important in rather old metal-poor stars (see e.g. 
Profhtt & Vandenberg 1991). We calculated models in- 
cluding microscopic diffusion and checked that in the 
Hyades, which are rather young, metal-rich objects, 
the effects of microscopic diffusion are small enough to 
be neglected in the present study. 

7. External boundary conditions. The r(T)-law is taken 
from the classical Eddington grey model atmosphere 
except in a few models where we use d, for compari- 
son, a r(T)-law derived from Kurucz's ( 1991 ) complete 
ATLAS9 model atmospheres. 

To fit the Hyades empirical M-L relation and HR di- 
agram, we calculated evolution models and isochrones 
with various values of the initial helium abundance Y 
and metallicity Z/X. The Z/X-ratio can vary in the 
range 0.034 ± 0.005 given by the observed [Fe/H]-value 
(see Sect. ^). The error bar on Z/X (0.005) includes 
the 0.05 dex error on [Fe/H] and the 11 per cent er- 
ror on {Z/X)q quoted by Anders & Grevesse ( 19891 ). 
The helium mass fraction Y is allowed to vary in the 
range 0.235-0.30, i.e. from the primordial abundance (see 
Peimbert et al. 2000[ ) to the approximative upper limit 
given by observations (see Nissen 1976). In the following, 
we refer to solar-scaled helium when the helium abun- 
dance is related to metallicity through the well-known re- 
lation iY-Yp)/Z^AY/AZ =(Ar/AZ)0where rp=0.235 
is the primordial helium value and {AY/ AZ) qIS derived 
from the calibration in luminosity and radius of the so- 
lar model. For a solar model that does not take into ac- 
count the microscopic diffusion we obtain Yq — 0.2674 
and Zq = 0.0175 and in turn (Ay/AZ)Q=1.9. 

The conversion of the theoretical model outputs (Mboi, 
Teff) into the observational plane (My, B — V) has been 
performed with the most recent version of the BaseL 
Stellar Library (BaSeL, version 2.2), available electroni- 
cally at http : //www, astro .mat .uc .pt/BaSeL/^ (Lejeune, 
2001). BaSeL provides color-calibrated theoretical flux 
distributions for the largest possible range of fundamental 
stellar parameters. Toff (2000 to 50,000 K), log g (-1.0 
to 5.5 dex), and [Fe/H] (-5.0 to +1.0 dex). 

The BaSeL flux distributions are calibrated on the 
stellar UBVRIJHKL colors, using these empirical pho- 
tometric calibrations for solar metallicity, and semi- 
empirical relations constructed from the color differences 
predicted by the stellar model atmospheres for non-solar 
metallicities (details about the calibration procedure are 
given in Lejeune et al. 1997, 199S). 

We point out that although Hipparcos provided its own 
precise magnitudes Hp (see Sect. 2J.), it has not been 



5. Constraints from the M-L relation 

The vB22-components have the lowest and most accurate 
masses (ctm/M < 1%) and define the lower part of the M- 
L relation quite accurately. From the models we find that 
at 65 Myr (age of the Hyades, according to Ferryman 
et al. 1998), both components are only slightly evolved. 



the departure from the ZAMS position being less than 
0.06 mag for vB22A (0.03 mag for vB22B). 

The positions of the vB22 components in the M-L 
plane are plotted in Fig. ^ where we have superimposed 
models and isochrones aged 625 Myr. The first important 
conclusion that can be drawn is that (with the present set 
of models) the constraints imposed by the M-L relation 
are hardly compatible with a solar-scaled helium value: 
the ([Fe/H], Y) = (0.14, 0.28) isochronc lies more than 
3cr above the data. This possibility had already been sus- 
pected by Torres et al. (|1997b|). Even if we allow [Fe/H] to 



possible to use them here because there are still problems 
in computing the bolometric corrections corresponding to 
the Hp band (Cayrel et al. |1997| ). 

All the models and isochrones presented here can be 
sent on request. 



vary inside its error bars, we find that an isochrone with 
[Fe/H]=:0.10 dex and solar-scaled helium (not represented 
in the figure) sits more than 2fj above the data. Errors on 
the bolometric corrections (BC) could be invoked but they 
should be quite high. The BaSeL transformations applied 
to (Tcff, log g)-values close to those of the vB22 system 
lead to BCa « —0.1 mag and BCs ~ —0.7 mag. On the 
other hand, for both stars, the BaSeL BCs do not differ 
by more than 0.03 mag from Alonso et al. 's (1995) em- 
pirical corrections. To reconcile models with solar-scaled 
helium with observations, the BC should be changed by 
ABCa ~ 0.15 mag and ABCs ^ 0.40 mag, which ap- 
pears to be very high. A helium content lower than the 
solar-scaled value is therefore favoured. 

The well-known degeneracy between helium and 
metallicity in the HR and M-L diagrams makes isochrones 
with different ([Fe/H], Y) values coincide exactly. Models 
give dMv/dY = -10 ± 1 and dMv/d[Ye/B\ = 2.0 ± 0.2. 
That means that a change in Y of —0.01 is compensated 
for by an increase in [Fe/H] by 0.05 dex. 

Several error sources bring uncertainties in the Y de- 
termination. The errors on mass and visual magnitude 
each lead to AY « 0.005. An error of 0.05 mag on the bolo- 
metric correction also gives AY « 0.005. Depending on 
the choice made for the Hyades [Fe/H]-value (i.e. spectro- 
scopic or photometric determination) the [Fe/H] contribu- 
tion to the y-error budget will be dominant or not: spec- 
troscopy gives A[Fe/H]=0.05 dex and in turn AY = 0.01 
while photometry yields A[Fe/H]=0.013 dex leading to 
AY = 0.0026. 

We also investigated how the models change when the 
input physics or input parameters are modified: 

— the component magnitudes are not changed by more 
than 0.001 mag if microscopic diffusion or pre-main 
sequence evolution calculations are included, 

— adopting the OPAL EOS instead of CEFF decreases 
both magnitudes by 0.02 mag, 

— changing the Eddington's grey model atmosphere into 
Kurucz's model only modifies the vB22B magnitude 
(by +0m mag). 
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— small age variations do not change the magnitudes 
much [dMv/dt < 10^^ mag per Myr), 

— varying ckmlt around Q!mlt,o yields dMy/daMLT — 
0.13 for vB22A (0.20 for vB22B), 

Using the OPAL EOS, we varied Y and omlt to find 
models of vB22A and vB22B satisfying = (A/y. model — 
Mv,obs)V'^Mv^ + (Mmodci - Mobs? /au^ < 1- We took ac- 
count of an error on the bolometric corrections of ABC — 
0.05 mag. Figure |^ shows the x^-contour levels in the Y- 
ctMLT plane. If we force auur to be in the range 1.2-2.0 
(i.e. Q!MLT=aMLT,©± 0.4, See item 4, Sect. ||) and if we re- 
quire that the initial helium content be the same in the two 
stars, we find solutions for Y in the range 0.247-0.263. In 
that y-range we can have auur, A — «mlt,b but the small- 
est are obtained for Aomlt = Q;mlt,a — Q!mlt,b 7^ 
with a trend for Aq;mlt to increase when Y decreases. 

In summary, the helium content of the Hyades deduced 
from models based on the OPAL EOS is F = 0.255±0.009 
([Fe/H] from photometry) or Y = 0.255 ± 0.013 ([Fe/H] 
from spectroscopy), that we write 

Y = [0.255 + 0.2([Fe/H] - 0.14)] ± 0.008, 

to highlight the effect of the [Fe/Hj-uncertainty. This 
value is fully compatible with Y= 0.26 ± 0.02 derived by 
Ferryman et al. ( |1998D from the model fit of 40 ZAMS 
stars with spectroscopic Toff. The changes in the input 
physics considered here do not modify that result by more 
than Ay = 0.002. Present models also suggest that omlt 
could be different in the two components, and higher in 
vB22A, but the error bars on observations are still too 
large to draw firm conclusions on that poin t. 

As pointed out by Lastennet et al. ( 1999| ), models with 
solar-scaled helium and a solar mixing-length ratio can- 
not satisfy the mass-radius (M-R) and the M-L relation 
simultaneously. From our models we find that, quite in- 
dependent of the helium abundance, only vB22A models 
with high ttMLT (Q!mlt>1-8) and vB22B models with very 
low ttMLT (aMLT<l-0) can give the observed radii. The si- 
multaneous fit of the M-L and M-R relations is possible 
for vB22A but imposes hardly acceptable values of q;mlt 
in vB22B. The use of Kurucz's model atmospheres to fix 
the external boundary conditions of the models does not 
modify this conclusion. 

New Toff-determinations by Cayrel de Strobel et al. 
( |200lD give Toff.A = 5600 ± 75 K and Teff,B = 
4500 ± 75 K. The radii of vB22A and vB22B can 
be derived from the observed magnitudes and these 
Teff through Stefan-Boltzmann's law and using bolomet- 
ric corrections from BaSeL. We find Ra — 0.97 ± 0.05 Kq 
and Rb = 0.64 ± 0.06 Rq. The radii based on Stefan- 
Boltzmann's law are hardly compatible with and much 
less precise than the radii inferred from observed eclipses 
(see Table 0). A higher Ra (smaller Rb) implies a smaller 
omlt.a (higher aMLT.e) and therefore it is easier to find 
models with Q!mlt in the range aMLT,0±O.4 that fit both 
the M-L relation and radii from Stefan-Boltzmann's law. 
To draw definite conclusions on the validity of the models. 



it is important to rediscuss the observations thoroughly. 
In particular, the radii from eclipses should be firmly as- 
sessed because they are given with such high precision that 
they strongly constrain the models. 

Figure |^ shows the whole M-L relation for the 5 bi- 
nary systems. Because the error bars on the component 
masses are quite large (i.e. in the range 5-25 percents), 
no additional constraint on the helium abundance comes 
from the four systems (/)342, 51 Tau, 9^ and 9^ Tau. In 
particular, 6*^ Tau A lies in the hook region (end of the 
core H-burning phase) but its exact position on the hook 
cannot be determined because of its large mass error bar. 

6. Revisited analysis of the HR diagram 

Figure || shows the colour- magnitude diagram of 90 secure 
Hyades members (de Bruijne et al. 200C| ) complemented by 
the data of Table ^ for four of the binary systems examined 
in Sect. || above {0^ Tau has been omitted because no 
reliable information on the individual {B — y)-values is 
available). A model isochrone, aged 650 Myr, calculated 
with the metallicity of the Hyades ([Fe/H]=0.14 dex) and 
a helium abundance fixed by the M-L relation (y=0.26) 
is also plotted in Fig. ^ showing that: 

1. In the range [B - V) G [0.3, 1.15] the fit of the data 
is quite good, although for {B — V) S [0.7, 0.9] some 
stars sit above the isochrone. 

2. For {B ~ V) > 1.15 the isochrone is much too blue 
with respect to the data: there is a change of slope at 
low mass both in the data and in the isochrone but it 
occurs at lower {B — V) in the isochrone. 

3. For {B — V)<0.3, corresponding to the turn-off region, 
a few stars are located to the right of the isochrone. 

Figure ^ shows the same data with, in addition, 54 
known (unresolved) binaries with dynamical parallaxes 
from Dravins et al. (1997). It shows that, except at low 
mass {{B — V) > 1.15) where the slope of the isochrone is 
bad, the binaries sit well in the AAfy ^ 0.75 mag band, 
where they are expected to be. 

In the following, we shall only discuss the MS and turn- 
off regions. The isochrones presented here do not go far 
enough to cover the giants, in the He-burning stage. In this 
region, only two secure members are found, plus 9^ Tau 
A. For a discussion of the constrain ts ass ociated with the 
giant region , see de Bruijne et al. (200C) and Castellani 
et al. (I2OOOI) . 



6.1. Helium abundance 

Figure ^ illustrates the differences between two isochrones 
of same age and same metallicity ([Fe/H]=0.14) but 
different helium abundances (F=0.26 and solar-scaled 
y=0.28). The isochrone that gives a better fit to the M-L 
relation (i.e. y=0.26) also provides a better fit to the low 
MS. However the constraints in the (My, B — V) plane are 
not so strong (the difference in B — V between the two 
isochrones are lower than the mean error bar on B — V) . 
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Like Castellani et al. (2000), we find that the isochrone 
with solar scaled helium departs from the data as soon as 
(B-V) > 0.6. Castellani et al. suggest that a decrease of 
(Xmlt could improve the fit at least for {B — V) £ [0.7, 1.1]. 
According to their Fig. 2 a rather large decrease of omlt 
would be necessary (by about ^^0.4). We discuss Q!mlt 



changes in Sect. 6.3 



6.2. Metallicity 

Fig. |8| shows two isochrones with extreme [Fe/H]-values 
(0.09 and 0.19) that were indistinguishable in the M-L 
plane (Sect. In the colour magnitude diagram they dif- 
fer, but by a very small amount. On the MS, a variation 
of 0.05 dex in [Fe/H] changes {B - V) by less than 0.01 
mag. The present data have {(T(^b-v)) — 0-01 mag which 
hardly allows to discriminate between [Fe/H]-values in- 
side the observed range, A[Fe/H]=0.013 dex (0.05) from 
photometry (spectroscopy). On the other hand, an uncer- 
tainty of 0.05 dex on [Fe/H] gives a 25 Myr uncertainty 
on the age. 

6.3. MS slope 

The effects of a change in auur on the (B—V) colour index 
are plotted in Fig. ^as a function of Aly for different values 
of mctalhcity (metal-poor [Fc/H]=:-1.0, solar [Fe/H]=0.0 
and Hyades [Fe/H]=+0.14 dex). We considered a ±0.20 
deviation of q;mlt from the solar value (q!mlt,0~ 1.6). At 
any metallicity a significant shift in {B — V) is expected 
for magnitudes in the range 4.-7. mag, with a maximum 
shift at My— 5.0-5.5 mag. The higher the metallicity, the 
higher the (B-V) shift. At Hyades metaUicity, the (B-V) 
shift is in the range ~0. 01-0. 03 mag. This effect is of the 
order of, and even larger than, the mean error bar on 
{B — V) for the Hyades data {{a(^B-v)) — 0-01 mag). 

Fig. ^is a zoom of the Hyades MS over the magnitude 
range where the omlt changes have the largest effects. 
Four isochrones with q;mlt=Q!mlt,0— 0.40, qimlt,©— 0.20, 
c>!MLT,0, and q;mlt,0+O.2O are superimposed on the data. 
All stars with {B - V)e [0.4, 1.1] he in the band where 
Qmlt, 0—0.4 <Q;MLT<aMLT,0+O.2. With the present set 
of model isochrones, we note a trend for stars of decreasing 
mass to lie on isochrones having decreasing ofMLT-valucs. 
This suggests aMLT-values in the range ^1.6-1.8 for 
M - 1.3 - 1.7 Mq, aMLTwl-4-1.6 for M - 1.0 - 1.3 Mq 
and ttMLT^ 1-4 below 1.0 Mq. 

6.4. Lower part of the MS 

There is a poor fit of the observed data in the cooler re- 
gions of the HR diagram (i.e. for {B — V) >1.2 mag, T'eff< 
4300 K). The discrepancy in {B — V) between models and 
observations amounts to about 0.07 mag (0.10 mag) at 
{B - y)-1.3 mag (1.4 mag). 

In the low-mass MS star region both the observations 
and the models presented here become less secure. Errors 



in the equation of state or external boundary conditions 
may have an important impact on the models. On the 
other hand, the B — V color index is not the best index 
to be used in the low mass red dwarf region. However, we 
found the same kind of discrepancy in other open clus- 
ters both in {B - V) and (V - /) (Robichon et al. |1999 , 
Lebreton pOOO| ). We now examine the input physics of the 
models. 



6.4.1. Equation of state 

In most models, we used the CEFF EOS (see Sect. |). 
Sophisticated EOSs, including collective effects, have been 
desig ned to study the Sun (MHDEOS, Mihalas et al. 
1988| ; OPAL EOS, Rogers et al. |1996D and very-low mass 



stars and planets (as the SC EOS, Saumon & Chabrier 
1991). The effects of the EOS on the models are illustrated 



in Fig. ^ where we compare models calculated with the 
CEFF and OPAL EOSs. Below 0.65 Mg, the MS slope is 
steeper in models including the most sophisticated EOSs 
and the fit with the observed low MS position becomes 
worse. We also calculated models with the MHD and SC 
EOS which exhibit the same behaviour. 



Lebreton & Dappen ( 1988 ) have shown that the M- 
L relationship is very sensitive to the EOS used at low 
mass. We calculated models with CEFF, MHD, OPAL 
and SC and found that at solar metallicity a model of 
~0.50 Mq is «0.75 mag brighter when calculated with the 
MHD, OPAL or SC EOS instead of the CEFF EOS. In the 
Hyades, no binaries have yet been observed in this mass 
range, but in the future new constraints for the physics 
will certainly come from low-mass binaries. 

6.4.2. Model atmospheres 

Model atmospheres are fundamental because they fix the 
external boundary conditions of the interior models and 
intervene in the transformations of the model results 
from the (Mboi, Toff)-plane to the {My, B — V) colour- 
magnitude diagram. 

The BaSeL conversions we used to convert the (Mboi, 
Tcff ) values to (My, B-V) are based on Kurucz's ATLAS9 
model atmospheres down to Toff~3500 K and on Allard 
& Hauschildt (|1995| ) non-grey model atmospheres beyond. 
Kurucz's ATLAS9 models are known to be insufficient be- 
low 4500 K but from preliminary tests we do not expect 
that conversions based on the new model atmospheres of 
Hauschildt et al. ( |l999[ ) will improve the fit in the 4000- 
4500 K range. Castellani et al. ( 2000| ) examined three sets 
of conversions (either empirical or based on Kurucz's mod- 
els) and did not find improvement of the slope-fitting at 
high {B — V). We also examined various empirical conver- 
sions and reach the same conclusions. 

On the other hand, the boundary conditions of 
our models are derived from the T(r)-laws based on 
Eddington's grey m odel a tmosphere. As demonstrated by 
Chabrier & Baraffe ( |l995 ), because convection penetrates 
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into optically thin layers in the envelope of low-mass stars, 
it is not correct to use grey model atmospheres and related 
T(T)-la"ws to fix the boundary conditions of the interior 
models for reff< 4000 K. Chabrier & Baraffe find that, at 
solar metallicity and in the range ~3000-4000 K, models 
based on the Eddington T(r)-laws are hotter and brighter 
than models based on non-grey full model atmospheres, 
the Teff-difference is of ~50 K at Tcff~3800 K and reaches 
~200 K at Toff~ 3300 K (then figure 2). 

We have here a discrepancy in {B — V) that amounts 
to 0.07-0.10 mag at B - y~1.3 mag (1.4 mag). From the 
BaSeL transformations we estimate that at Mv~8., the 
effective temperature of our models (which is ~4150 K) 
should be reduced by more than 250 K to increase the 
(B — V) index to the observed position. This appears to 
be very large with respect to Chabrier & Baraffe's results. 

At high temperatures (>4500K), models calculated 
with outer boundary conditions taken from T(r)-laws 
based on Kurucz's model atmospheres and models us- 
ing Eddington's grey law differ by less than 0.01 mag in 
(B-V). 

6.5. Turn-off region and age 

The turn-off region in the Hyades corresponds to Afy- 
values below 2.5 mag. As shown in Fig. ^ there are only 
12 turn-off stars among the secure stars selected by de 



Bruijne et al. (2000). We have added to the figure the two 
components of 9'^ Tau and the brighter component of the 
51 Tau system. 

The interpretation of observations in the turn-off re- 
gion of the Hyades are complicated by the effects of rota- 
tion and overshooting that make either models or photo- 
metric data uncertain. 

6.5.1. Overshooting 

Hyades stars at turn-off have masses in the range 2.0- 
2.5 Mq. Therefore, they have had convective cores during 
their MS phase. Overshooting of convective cores extends 
the mixing in central regions: it increases the amount of 
hydrogen available for H-burning and in turn, increases 
the MS lifetime and the MS width in the HR diagram. 
The amount of overshooting is poorly known. Recently, 
Ribas et al. ( 2000D proposed a mass dependence of the 
overshooting on the basis of the modeling of 8 eclipsing 
binaries. At M ~ 2. M© and for solar metallicity, they 
give dov — 0.17±0.05Hp which is fully c ompat ible with 

dov -- 
here 



0.2Hp, the value of Schaller et al. ( |1992| ) adopted 



As discussed in Ferryman et al. (1998) who compared 
two sets of isochrones with and without overshooting, it 
is possible to fit the turn-off in both cases at different 
ages. This is shown again in Fig. 12 where we have plotted 
isochrones of various ages with and without overshooting. 
The theoretical isochrones without overshooting are the 
same as in Ferryman et al. (1998) but converted into the 



(My, Toff) plane by means of the BaSeL library. For the 
same [B — V] at turn-off, models without overshooting are 
about 100 Myr ^ younger than models with dov = 0.2Hp. 
Models without overshooting appear to give a poorer fit of 
the position of fl^Tau A, and at the magnitude of ^^Tau A, 
they predict a mass 2 percent larger than the mass pre- 
dicted by models with overshooting. However, the mass 
of 0'^ Tau A is not known with enough precision to dis- 
criminate between models with dov = 0.2Hp and models 
without overshooting. Furthermore, as will be discussed 



in Sect. 6.5.2, one must also account for rotation effects in 



the Hyades turn-off region. 
6.5.2. Rotation 

The turn-off region of the Hyades also corresponds to the 
5 Scuti instability strip. In that region, stars are mainly 
rapid rotators and variable stars and, some of them have 
been identified as (5-Scuti pulsators (Antonello & Pasinetti 
Fracassini 19981 ). The 15 stars seen in Fig. have Ve sinz 
in the range 4 3-175 km.s^i (Bernacca & Perinotto |1970| , 
Abt & MorrcU 1995) and among these stars, 10 stars have 
Vf, sm I > 90 km.s^i (ue is the equatorial velocity and i 
the inclination of the rotation axis with respect to the 
line of sight). In particular, 51 Tau A has measured values 
of Ve sini in the range 97-125 km.s~^, 6^ Tau A (which is 
a well-known (5-Scuti) has UgSini ^ 80 km.s^^ while its 
companion 9^ Tau B rotates even faster with published 
WeSini in the range 90-170 km.s~^ (Torres et al. 1997c ). 

Rotation complicates the interpretation of the HR di- 
agram and makes the age estimates uncertain. It affects 
both the models and the photometric data. 



It has been known for a long time that rotation pro- 
duces a displacement of the stars in any photomet- 
ric HR diagram by amounts depending on and i 
(Maeder & Peytremann 1972). Perez-Hernandez et al. 
(1999) have computed the photometric magnitudes of 
rotating atmospheres. They give the magnitude differ- 
ences between non-rotating stars of given (Tcff, log g) 
and their rotating copartners, as a function of the i- 
angle and oi lo = Q./Q,c where f2 is the star angular 
velocity and = (8GM/27i?3)i/2 {^m is the mass, 
Rp is the polar radius). Interpolation in their Table 2 
at (Toff, log g)-values of the model isochrones yields 
the corrections AMy and A{B — V) to apply to 
our My and [B — F)-values. We calculated the cor- 
rections for three values of WeSini (50, 100 and 150 
km.s~^), and for each iigSini, we considered the case 
where i=90 degrees {ve and therefore w are minimum), 
and the case where ld = 0.95 (ug tends to the maximum 
allowed value, i tends to its minimum possible value). 



Our models are based on ~700 mass shells. Breger 



Pamyatnykh (199J) pointed out that the MS width is reduced 



if models are calculated with smaller time and space steps, the 
effect being larger in models with overshooting. This numerical 
effect could result in an overestimate of the ages in the turn-off 
region of ~50 Myr. 
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The results are plotted in Fig. |l3|. When the rotation 
rate is close to maximum, the isochrone is shifted both 
toward higher luminosities and smaller effective tem- 
peratures. The shift in {B — V) is of a few tenths of a 
magnitude and does not depend much on Ve sini, ex- 
cept in the upper part of the turn-off. For minimum ro- 
tation rates (i=90 degrees), the effect of rotation is to 
shift the isochrones, mainly in (B — V), by an amount 
that increases with Vg . Rotating stars are therefore ex- 
pected to be spread in a band that extends on the 
red side of the non-rotating isochrone. A detailed in- 
dividual analysis of each star at turn-off is beyond the 
scope of this work. However, with the help of Fig. |l^, 
we can estimate that the effect of rotation on photom- 
etry could be responsible for an overestimate of the 
age of the Hyades by ~50 Myr. 

On the other hand, models are modified when rotation 



is taken into account. Breger & Pamyatnykh ( 1998 ) 
showed that models with uniform rotation are dis- 
placed toward the right in the HR diagram, in the same 
way as models including overshooting. Furthermore, 
the rotation profile inside a star is related to the re- 
distribution of the angular momentum by internal mo- 
tions which could be generated by meridional circula- 
tion and shear turbulence in a rotating medium (see 
Zahn 1992). Such motions might induce internal mix- 
ing. As shown by Talon et al. (1997), in the HR dia- 
gram these rotational effects also "mimic" overshoot- 
ing (for instance in a 9 Mq star with a rotational 
velocity of 100 km.s^^ the effect is roughly equivalent 
to the effect of an overshooting of 0.2Hp). 



Overshooting and rotation have similar effects on 
model isochrones and therefore similar signatures in the 
HR diagram. These effects cannot be discriminated with 
the present data; additional data on the internal struc- 
ture of the considered stars are required that will hopefully 
(and only) be brought by asteroseismology measurements. 
On the basis of the present data and model isochrones, the 
upper limit on the Hyades age is around ~ 650 Myr, the 
lower limit would be ~500 Myr (no overshooting, rotation 
corrections on photometry included). 

7. Summary and conclusions 

High accuracy observations including astrometric, photo- 
metric and spectroscopic data have been gathered for the 
nearest open cluster, the Hyades. These data provide pre- 
cious information: the Hyades is the only cluster in which 
individual distances are available for a bunch of a hundred 
stars and in which individual masses have been measured 
in several binary systems. These strong constraints have 
been used to infer the characteristics of the cluster from 
the comparison with stellar models. 

The helium abundance of the cluster can be con- 
strained by the M-L relation provided (1) the stellar data 
are accurate enough, (2) the error on the bolometric cor- 
rections needed to convert the model results in the ob- 



served M-L plane are small and, (3) the input physics 
of the models is well chosen. The only Hyades binary 
system with data accurate enough to constrain the he- 
lium abundance is vB22. We estimated the vB22 helium 
abundance using models including different input physics 
(EOS, atmospheres, diffusion, PMS) and we allowed vari- 
ations of ttMLT of ±0.4 around the solar value. The result, 
Y — 0.255 ± 0.009, implies that the Hyades have less he- 
lium than expected from their supersolar metallicity. A 
low helium abundance has already been suggested in the 
past ( see fo r instance VandenBerg & Poll 



198C, Torres 



et al. |l997q) and, a long time ago, on the basis of pho- 
tometric measurements, Hyades stars were found to be 
helium deficient with respect to field stars of same metal- 
licity (Stromgren et al. 1982). 

We pointed out that the error bar on Y is much smaller 
if [Fe/H] is taken from photometr y (int ernal error bar on 
[Fe/H] around 10 percent, Grenon ^OOOj ) rather than from 
spectroscopy (internal error bar on [Fe/H] around 35 per- 
cent, Cayrel et al. 1997| ). Progress will come from a better 
description of the so-called metallicity by the future deter- 
minations of the individual abundances of the main metal- 
licity contributors. The C, N, O abundances are presently 
poorly known while the Ne abundance is unknown (R. 
Cayrel, private communication). On the other hand, it 
would be worthwhile to improve the mass values in the 
four other binary systems which are presently not good 
enough to constrain models (and helium) significantly. 

An upper limit of ^^650 Myr on the cluster age can be 
estimated from the positions of stars at turn-off. It is diffi- 
cult to estimate precisely the age of the Hyades because at 
turn-off rapidly rotating stars are found which had a con- 
vective core on the MS. Rotation affects the photometry 
and makes the observed HR diagram positions uncertain, 
whereas internal rotation and overshooting are responsi- 
ble for internal mixing that may change the model course 
in the HR diagram. As a result, the age of the Hyades 
might be in the range 500-650 Myr. An improvement of 
the mass of 0^Tau A would certainly better constrain the 
overshooting by anchoring the star more precisely on its 
isochrone. Asteroseismological measurements and analysis 
are required to get a better understanding of the transport 
processes inside stars (Goupil et al. 200C| ). 

The observed MS slope, in the region of the HR di- 
agram corresponding to masses in the range 0.8-1.5 Mq, 
is quite well-defined and provides first tests of the auur- 
value to be used in models. It suggests that auLT could be 
a slowly decreasing function of mass: umlt would decrease 
from values of ~ 1.8- 1.9 at high mass to values of '--^ 1.2- 1.4 
at low mass. The confirmation of that trend would require 
extremely precise radii or Tes for individual stars along the 
MS. Presently, both radii and T^s are only available for 
vB22 (only vB22A sits in the HR diagram region sensitive 
to ckmlt) and since they do not agree, we have not been 
able to draw firm constraints on omlt from these data. 

In the very low mass region of the HR diagram, models 
are definitely too blue. This region provides tests on the 
model external boundary conditions and on the equation 
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of state and the knowledge of a few masses would help to 
understand the origin of the discrepancy. 

The Hyades stars are now better understood because 
they are fully seen as individuals and because of the 
considerable recent progress in the description of the 
physics of the stellar plasma (opacities, EOS). In 1989, 
VandenBerg & Poll (198!;) used stellar models to infer the 
Hyades distance modulus by means of the MS fitting tech- 
nique assuming that the helium content was solar. Now 
that we have removed the distance problem, we are able 
to test finer details in the models. We can hope that in 
the decade to come, a substantial number of open clusters 
will come to be known at the level of accuracy presently 
reached for the Hyades and even better. Hipparcos pro- 
vided mean distances to a few open clusters which already 
give rise to many questions. However, the discussion of the 
helium abundance in the nearer clusters is limited because 
errors sources can be important (color-calibrations and 
bolometric corrections, distances, metallicity) and because 
of the lac k of w ell observed binaries (Robichon et al. 1999, 
Lebreton 2000 ). In the future, GAIA will determine indi- 
vidual distances of stars in many clusters together with 
information on their abundances (Ferryman et al. 1997) : 
it is expected that about 120 clusters, located to within 
1 kpc, will reach to the level of accuracy or even better 
than what is reached today in the Hyades. GAIA will also 
discover binary systems in which it will measure masses 
accurately (at least ten binary systems are expected per 
cluster). Other observational devices will also provide ac- 
curate masses and radii in the near future. Therefore, we 
can expect that we shall go on constraining the models to 
progress in the understanding of the chemical and dynam- 
ical evolution of our Galaxy. 

Acknowledgements. We warmly thank D. Dravins, L. 
Lindegren and J. de Bruijne who provided their parallax 
data. We are grateful to A. Baglin, G. Cayrel de Strobel, 
R. Cayrel, H.-G. Ludwig, F. Arenou, M.J. Goupil and E. 
Michel for many fruitful discussions. We thank the referee Dr. 
A. A. Pamyatnykh for interesting remarks and information. 
Y. Lebreton acknowledges University of Rennes 1 for working 
facilities. This work was partially supported by the "Convnio 
ICCTI Embaixada de Frana, number 060/BO" and by the 
project "ESO/P/PRO/12128/1999" from "Fundao para a 
Cincia e Tecnologia". 

References 

Abt, A.H., MorreU, N.I., 1995, ApJS 99, 135 
Alexander, D.R., Ferguson J.W. 1994, ApJ 437, 879 
AUard, F., Hauschildt, 1995, ApJ 445, 433 
Hauschildt, P. H., AUard, F., Baron, E., 1999, ApJ 512, 377 
Alonso, A., Arribas, S., Martmez-Roger, C. 1995, A&A 297, 
197 

Anders, E., Grevesse, N., 1989, Geochim. Cosmochim. Acta, 
53, 197 

Andersen, J. 1991, AAR 3, 91 

Antonello, E., Pasinetti Fracassini, L.E., 1998, A&A 331, 995 
Bernacca, P.L., Perinotto M., 1970, A catalogue of stellar ro- 
tational velocities, Contr. Oss. Astrof. Padova in Asiago, 
239 



BesseU, M.S., Castelh, F., Plez, B. 1998., A&A 333, 231 
Bohm-Vitense, E. 1958, Z. Astrophys. 54, 114 
Breger, M., Pamyatnykh, A.A. 1998, A&A 332, 958 
de Bruijne, J.H.J., Hoogerwerf, R., de Zeeuw, P.T. 2001, A&A 
367, 111 

Castellani, V., degl'Innocenti, S., Prada Moroni, P.G. 2001, 

MNRAS 320, 66 
Caughlan, G.R., Fowler, W.A. 1988, Atomic Data Nuc. Data 

Tables 40, 284 

Cayrel, R., Castelli, F., Katz, D., van't Veer, C, Gomez A.E., 
Perrin, M.-N., 1997, in Hipparcos Venice '97, ESA SP-402, 
pp. 433 

Cayrel de Strobel, G., Crifo, F., Lebreton, Y., 1997, in 

Hipparcos Venice '97, ESA SP-402, pp. 687 
Cayrel de Strobel, G., Cayrel, R., Lebreton, Y., Morel, P., 2001, 

A&A in preparation 
Chabrier, G., Baraffe, L, 1995, ApJ 451, L29 
Christensen-Dalsgaard, J. 1991, in 'Challenges to theories of 

the structure of moderate-mass stars'. Cough D.O. & 

Toomre J. (eds.). Springer- Verlag, 11 
Dravins, D., Lindegren, L., Madsen, S., Holmberg, J., 1997, in 

Hipparcos Venice '97, ESA SP-402, 733 
Eggleton P.P., Faulkner J., Flannery B.P. 1973, A&A 23, 325 
Fernandes, J., Lebreton, Y., Baglin, A. et al. 1998 A&AS 338, 

455 

Goupil, M.-J., Dziembowski, W. A., Pamyatnykh, A. A., Talon, 
S., 2000, in "Delta Scuti and Related Stars", Eds. M; 
Breger & M. Montgomery, ASP Conf. Ser., Vol. 210. pp 
267 

Grevesse, N., Noels, A., 1993, in 'Origin and Evolution of the 
Elements', Prantzos N., Vangioni-Flam E., Casse M. (eds.), 
Cambridge: Cambridge University Press 

Grenon, M., 2000, in 'Luminosity Calibrations with Hipparcos', 
Highlights of Astronomy Vol. 12, Joint Discussion 13 of the 
XXIIIIrd General Assembly of the lAU, T.L. Evans (ed.), 
in press 

Iglesias, C.A., Rogers, F.J. 1996, ApJ 464, 943 

Kurucz, R.L. 1991, in 'Stellar Atmospheres: Beyond Classical 

Models', ed. L. Crivellari, I. Hubeny, D.G. Hummer, Nato 

ASI Ser., pp. 441, Dordrecht: Kluwer 
Lastennet, E., Valls-Gabaud, D., Lejeune, Th., Oblak, E., 1999, 

A&A 349, 485 

Lebreton Y., 2000, in Hipparcos and the Luminosity 
Calibration of the Nearer Stars, 24th lAU GA, JD 13, 
Highlights of Astronomy Vol. 12, T.L. Evans ed., Kluwer 

Lebreton, Y., Dappen, W. 1988, in Proc. Symp. Seismology of 
the Sun and Sun-like Stars, Tenerife, Spain, ESA SP-286, 
661 

Lebreton Y., Gomez A.E., Mermilliod J.-C, Ferryman M.A.C. 

1997, in Hipparcos Venice '97, ESA SP-402, p. 231 
Lebreton, Y., Perrin, M.-N., Cayrel, R., Baglin, A., Fernandes, 

J., 1999, A&A 350, 587 
Lejeune, T., 2001, in preparation 

Lejeune, T., Buser, R., 1999, in "Spectrophotometric dating 
of stars and galaxies", eds. I. Hubeny, S. Heap, and R., 
Cornett, ASP Conf. Series, Vol. 192, 211 
Lejeune, T., Cuisinier, F., Buser, R. 1997, A&AS, 125, 246 
Lejeune, T., Cuisinier, F., Buser, R. 1998, A&AS, 130, 75 
Ludwig, H.G., Freytag, B., Steffen, M. 1999, A&A 346,111 
Meader, A., Peytremann, E., 1972. A&A, 21, 279 
Martin, C, Mignard, F., Hartkopf, W.L, McAlister, H.A., 

1998, A&AS 133, 149 
McClure, R.D., 1982, ApJ 254,606 



Y. Lebreton et al.: The helium content and age of the Hyades 



11 



MermiUiod, J.-C, 1998, in 'Dynamical Studies of Star Clusters 
and Galaxies', Joint European and National Astronomy 
Meeting for 1998 (JENAM98), Kroupa P., Palous J. & 
Spurzem R. (eds.), Abstract Book published by ESA, The 
Netherlands. 

Mihalas, D., Dappen, W., Hummer, D.G., 1988, ApJ 331, 815 

Morel, P.J. 1997, A&AS 124, 597 

Narayanan V. K., Gould A., 1999, ApJ 515,256 

Nissen, P.E., 1976. A&A 50, 343 

Peimbert, M., Peimbert, A., Ruiz, M.T. 2000, ApJ 541, 688 
Perez-Hernandez, F., Claret, A., Hernandez, M.M., Michel, E., 

1999, A&A, 346,586 
Perryman, M., Brown, A.G.A., Lebreton, Y., Gomez, A.E., 

Turon, C, Cayrel de Strobel, G., MermiUiod, J.-C, 

Robichon, N., Kovalesky, J., Crifo, F., 1998, A&A 331, 81 
Perryman, M.A.C., Lindegren, L., Turon, C, 1997, in 

Hipparcos Venice '97, ESA SP-402, 743 
Peterson, D.M., Solensky, R., 1988, ApJ 333, 256 
Peterson, D.M., Stefanik, R.P., Latham, D.W., 1993, AJ 105, 

2260 

Proflitt, C.R., VandenBerg, D.A., 1991, ApJS 77, 473 
Ribas, L, Jordi, C, Gimnez, A., 2000, MNRAS 318, L55 
Robichon, N., Arenou, F., Lebreton, Y., Turon, C, MermiUiod, 

J-C, 1999, in Harmonizing Cosmic Distance Scales in a 

Post-HIPPARCOS Era, D. Egret & A. Heck eds., ASP 

Conference Series 167, pp. 72 
Rogers, F.J., Swenson, F.J., Iglesias, C.A., 1996, ApJ 456, 902 
Saumon, D., Chabrier, G., 1991 Phys. Rev. A 44, 5122 
Schaller, G., Schaerer, D., Meynet, G., Maeder, A., 1992, 

A&AS 96, 269 
Schiller, S.J., Milone, E.F., 1987, AJ 93, 1471 
Schwan, H., 1991, A&A 243, 386 
Sderhjelm, S., 1999, A&A 341, 121 

Strmgren, B., Olsen, E.H., Gustafsson, B., 1982, PASP 94, 5 
Talon, S., Zahn, J.-P., Maeder, A., Meynet, G., 1997, A&A, 
322, 209 

Torres, G., Stefanik, R.P., Latham, D.W., 1997a, ApJ 474, 256 
Torres, G., Stefanik, R.P., Latham, D.W., 1997b, ApJ 479, 268 
Torres, G., Stefanik, R.P., Latham, D.W., 1997c, ApJ 485, 167 
VandenBerg, D.A., PoU, H.E., 1989, AJ 98, 1451 
Zahn, J. P., 1992, A&A 265, 115 




1.09 



1.08 



1.07 



1.06 



1.05 




0.76 



Sun 



Fig. 2. The vB22 system in the M-L plane (data from 
Tables |l| & ||). The ellipses represent the locus where 
(Mv - Mv,obs)V<^Mv^ + (M - Mobs)Vf^M' < 1. Models 
aged 625 Myr are plotted. Continuous line: three indis- 
tinguishable isochrones with ([Fc/H], Y) = (0.09, 0.25), 
(0.14, 0.26) and (0.19, 0.27). Dot-dashed line: isochrone 
with [Fe/H]=0.14 dex and solar-scaled heUum y=0.28. 
Circles: models with ([Fe/H], Y) = (0.14, 0.26) and 
aMLT= aMLT,o-0-20 (•) or q;mlt.0+O.2O (o). Squares: 
([Fe/H], Y) = (0.09, 0.24), (0.14, 0.25), (0.19, 0.26) (□); 
([Fe/H], Y) EE (0.09, 0.26), (0.14, 0.27), (0.19, 0.28) (■). 
Open triangle: ([Fe/H], Y) = (0.14, 0.26) and OPAL EOS. 
The smallness of the error bars allows us to discriminate 
between rather close values of Y and to eliminate the 
isochrone with solar-scaled helium. 
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Fig. 3. Contour levels in the Y-auLT plane. Continuous 
lines are for = (My, model - My ,ohs)^ / cr Mv^ + 

(Mmodci - Mobs)V<7M^ = 10-2, dashed for = 1. vB22B 
is in bold lines. 




B-V 



Fig. 5 . Hya des 's HR diagram. Data are from de Bruijne 
et al. ( 2000| ) complemented by the data of Table || for four 
of the binary systems (o). Isochrone: 650 Myr, Y—Q.26, 
[Fe/H]=:0.14, aMLT=aMLT,0 



M„ 




2.0 1.5 
M/M, 

Sun 

Fig. 4. The Hyades M-L relation for vB22 (1), 51 Tau (2), 
0342 (3), 9^ Tau (4) and 9'^ Tau (5) and model isochrones 
aged 625 Myr. Continuous line: three indistinguishable 
isochrones with ([Fe/H], Y) = (0.09, 0.25), (0.14, 0.26) 
and (0.19, 0.27). Dot-dash line: [Fe/H]=0.14 dex and solar- 
scaled helium y=0.28. Here the error bars on mass are too 
large to discriminate between various helium values (vB22 
excepted) 



B-V 

Fig. 6. The HR diagram including unresolved known bi- 
naries. Open squares are data for single stars from de 
Bruijne et al. (2000), open circles are for the four resolved 
binary systems and full squares are the unresolved bina- 
ries with dynamical parallaxes from Dravins et al. ( 1997 ). 
Continuous line is the (650 Myr, r=0.26, [Fe/H]=0.14)- 
isochrone, dashed line is the same isochrone translated by 
—2.5 log 2 in My giving the maximum shift expected from 
unresolved binaries. The dot-dashed isochrone is 50 Myr 
old. 
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Fig. 7. Helium: zoom on the non-evolved region. Fig. 9. Effect of a ehange of ajviLT on the MS position. 
Isochrones aged 650 Myr, with [Fe/H]=0.14 and F=0.26 A{B — V) is the shift at constant magnitude between an 



(continuous line) or solar-scaled helium y=0.28 (dash) 




isochrone with aMLT=Q:MLT,0+AQ:MLT and an isochrone 
with aMLT,0- Positive values of A{B — V) correspond 
to Aq!mlt= —0.2, negative values to AaMLT= +0.2. 
The continuous line is for [Fe/H]=0.0, dot-dash line for 
[Fe/H]= -1.0 and dash line for [Fe/H]= +0.14 dex 
(Hyades) 



B-V 

Fig. 8. Two isochrones aged 650 Myr with "ex- 
treme" [Fe/H]-values are plotted: the continuous line 
is for ([Fe/H]=0.09, F=0.25), dashed for ([Fe/H]=0.19, 
F=0.27) 
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Fig. 10. Effect of a change of omlt on the MS po- 
sition at Hyades metahicity. Isochrones aged 625 Myr, 
[Fe/H]=0.14, Y=0.26, are plotted: the continous isochrone 
is for aMLT,0, dashed for aMLT^ckMLT,©— 0.2, dot- 
dashed for aMLT=ckMLT— 0.4, and dot-dot-dashed for 
aMLT=aMLT,0+O.2. The positions of models of different 
masses on the isochrones are indicated by arrows. Open 
circles are vB22A and vB22B 
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Fig. 12. Zoom on the turn-off region. The observation- 
nal data are the same as in Fig ^. Isochrones with 
([Fe/H]=0.14, y=0.26) and different ages are plotted. 
From left to right: continous lines are isochrones with over- 
shooting and aged 600 Myr, 625 Myr (bold), 700 Myr and 
750 Myr, dashed lines are isochrones without overshooting 
and aged 550 Myr (bold) and 600 Myr 
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Fig. 11. Effect of the EOS on the lower MS position. The 
continuous line is the isochrone calculated with the CEFF 
EOS; dashed line with the OPAL EOS. 
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Fig. 13. Effects of rotation on photometric data. 
Continuous lines are isochrones with from left to right 
Ve sini = 0, 50, 100, 150 km.s~^ and z=90 degrees. Dashed 
lines are ischrones with, from left to right (at Mv=0.5), 
VeSini = 50,100,150 krn.s"^ and oj = 0.95 (see text). 
Open (full) symbols are Hyades with Ve s'mi > 90 km.s^^ 
{ve sinz < 90 krn.s"^). Circles are for the binaries studied 
in Sect. ^ 



